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ABSTRACT 

The nonlinear stability of a compressible flow past a wedge is investigated in the hyper- 
sonic limit. The analysis follows the ideas of a weakly nonlinear approach as first detailed 
by Smith (1979). Interest is focussed on Tollmien-Schlichting waves governed by a triple 
deck structure and it is found that the attached shock can profoundly affect the stability 
characteristics of the flow. In particular, it is shown that nonlinearity tends to have a sta- 
bilising influence. The nonlinear evolution of the Tollmien-Schlichting mode is described in 
a number of asymptotic limits which were first identified by Cowley and Hall (1990) in their 
linearised account of the current problem. 
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§1. Introduction. 

This investigation is an extension of the work of Cowley & Hall (1990) who were 
concerned with the linear stability of hypersonic flow around a wedge of small angle. Here 
we aim to examine the effects of including nonlinearity in the solutions for the viscous 
modes of instability considered in Cowley Sz Hall. 

Recently, there has been a resurgence in research of hypersonic flows, motivated by 
the interest in developing high speed aircraft. In compressible flow work it is necessary to 
specify a functional relationship between the temperature and the viscosity of the fluid. 
Frequently Chapman’s law is used in which the viscosity is taken to be directly proportional 
to the temperature. Although this law can be rather unrealistic for high speed flows, 
stability analysis using this relationship is often more straightforward than for any other 
commonly used law. Cowley Sz Hall (1990), hereafter referred to as CH, used Chapman’s in 
their work and we shall follow their lead here. A much more accurate relation, Sutherland’s 
law, is used by Blackaby, Cowley & Hall (1990) in their investigation of the instability of 
hypersonic flow past a flat plate. They obtained smaller growth rates for the vorticity 
mode under consideration compared to those when Chapman’s law was assumed. In this 
instance we know that Chapman’s law is a bad approximation for the vorticity mode since it 
is concentrated in the layer where there is substantial variation in the basic temperature. 
However, since we are going to be concerned with the stability of Tollmien-Schlichting 
waves which by their nature are dominated by viscous effects close to the wall it is be to 
hoped that the results obtained here using Chapman’s law will not be vastly different to 
those found should a more accurate law be used. 

At hypersonic speeds real gas effects will be important but we do not take them into 
account in this study. These have been considered by Fu, Hall & Blackaby (1990) in 
their investigation of the instability of Gortler vortices at hypersonic speeds, in which they 
also assumed Sutherland’s law. The fluid was taken to be an ideal gas and the effects 
of gas dissociation on the stability of the modes considered were obtained. The Gortler 
instability at hypersonic speeds has also been investigated by Hall Sz Fu (1989) and Dando 
Sz Seddougui (1991). 

As already mentioned, in the present study attention is focussed on (viscous) Tollmien- 
Schlichting waves. CH investigated the effect of an attached shock on the viscous and 
inviscid modes in the compressible boundary layer on a thin wedge. For a discussion 
of the inviscid modes the reader is referred to Mack (1987). Smith & Brown (1990) 
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have investigated inviscid modes at high Mach numbers and in the absence of a shock. 
These authors concentrated on the vorticity mode while CH, in their study of the inviscid 
modes, looked at the acoustic modes which differ from the vorticity modes in as much as 
they are not concentrated within a thin layer. Blackaby, Cowley & Hall (1990) discuss 
the relationship which exists between these various modes. For large Mach numbers the 
inviscid modes have larger growth rates than the viscous ones but nevertheless the effect of 
the shock on the viscous modes is of some importance. In fact CH show that the presence 
of the shock greatly affects the viscous modes. 

The wavelengths and frequency of the viscous instability in a compressible boundary 
layer at hypersonic speeds were given by Smith (1989). These form the basis of the 
scalings adopted in the present analysis. CH described the interactive triple-deck structure 
governing disturbances to a hypersonic flow over a wedge. The basic flow is given by an 
exact steady solution of the governing equations for a constant density gas and we follow 
CH and assume for simplicity that the wedge is insulated. As mentioned by CH the 
analysis can be extended to consider an isothermal flow where heat transfer at the surface 
is allowed. CH were concerned with the stability of the basic flow to viscous modes and as 
a first analysis they sought conditions under which non-parallel effects could be neglected. 
They showed that in order to do this it is necessary to make the Newtonian approximation 
and this yields a lower bound on the Mach number. A discussion of when non-parallel 
effects become important for hypersonic flow is given by Smith (1989) and this gives the 
same lower bound for the Mach number. CH also selected the angle between the surface 
of the wedge and the shock so that the shock occurs in the upper deck of the triple-deck 
structure. In this way the influence of the shock wave on the Tollmien-Schlichting waves 
can be investigated. As in classical triple-deck analysis as the perturbation to the basic 
flow is increased in amplitude from an infinitesimal size to a point where nonlinearity must 
be accounted for, it is the equations governing the disturbance in the lower deck which 
are the first to become nonlinear. We choose to retain the nonlinearities in the lower deck 
equations since the aim of the present study is to investigate these effects. The above 
requirements give bounds on the position of the shock or alternatively bounds on the size 
of the local Mach number with respect to the size of the Reynolds number, which will be 
taken as large for the following analysis. 
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CH determined the linearised conditions which hold at the shock. Since the nonlinear 
effects are unimportant in the upper deck (where the shock occurs) to the order of approx- 
imation considered both by CH and the present study these linearised shock conditions 
are appropriate. It is found that the effect of the shock is weak. 

The linear results of CH for the Tollmien-Schlichting waves will be presented in §4, 
however, we give a brief summary of their results here. CH determined the nonlinear system 
governing the instability of Tollmien— Schlichting waves over a thin wedge in the presence of 
a shock. They linearised this system in order to determine the dispersion relation relevant 
to small disturbances and neutral solutions were presented for various positions of the 
shock. Since the disturbances are three-dimensional there exists an infinity of solutions. 
CH presented asymptotic solutions for the neutral wavenumbers which agree well with 
their numerical results. In particular, they showed that as the position of the shock, y s , 
tends to infinity the shockless dispersion relation derived by Smith (1989) is recovered. 
CH also gave a discussion on the temporal growth of the linearised modes. 

CH proceeded to investigate the effect of the shock on neutral, two-dimensional invis- 
cid modes of instability having wavelengths scaled on the boundary layer thickness. These 
modes are concentrated in the temperature adjustment layer centred at the generalised 
inflection point and solutions were given for the acoustic modes alluded to earlier. (A 
discussion of the vorticity modes in the absence of a shock is given by Smith & Brown 
(1990).) However, CH show that the shock has a negligible effect on the acoustic modes 
unless the distance of the shock from the surface is of the order of the boundary layer 
thickness. 

Before we give an outline of the paper we mention some additional recent investigations 
of instability of hypersonic flow. Balsa & Goldstein (1990) investigated the high Mach 
number instability of supersonic mixing layers using asymptotic theory. They use the 
WKB method and solve the compressible Rayleigh equation for a non-isothermal free 
shear flow, so in a sense their solutions are similar to CH. This work is also related to that 
on the vorticity mode by Smith & Brown (1990). An investigation of nonlinear effects 
on the acoustic modes at hypersonic speeds described by CH is given by Goldstein & 
Wundrow (1990). Many aspects of theoretical research into hypersonic flow are discussed 

in the recent paper by Brown et al. (1991). 

The plan of the paper is as follows. In §2 we describe the basic flow situation for 
hypersonic flow around a thin wedge, and give the equations governing the flow. The 
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scalings are chosen so that we are able to neglect non-parallelism, to include the effects of 
nonlinearity in the lower deck of the triple-deck structure and to ensure that the shock is 
located in the upper deck. The triple-deck structure for hypersonic flow is described in §3 
and the scalings chosen to eliminate some of the variables for convenience. In addition, 
we state the linear conditions which hold at the shock (and which are described in detail 
in CH). In §4 we consider the nonlinear solutions of the governing system of equations. 
Rather than attempt a full numerical solution of the nonlinear equations we execute a 
weakly nonlinear analysis of the system in a manner virtually identical to that of Smith 
(1979) who studied the nonlinear stability of an incompressible boundary layer over a 
flat plate. In this way we derive an amplitude equation of classical Stuart-Watson type to 
describe the evolution of the disturbance. Solutions of this evolution equation are discussed 
in §5 where we also draw some conclusions. 


§2. The basic flow. 

We consider the flow of a compressible fluid over a wedge of semi— angle 8. As in CH, 
we suppose that the supersonic flow has velocity magnitude U and we consider the situation 
when the wedge is symmetrically aligned with the oncoming flow. Then symmetrical shocks 
develop on either side of the wedge. We denote the semi-angle of the shocks by a and so 
the angle between the shock and the wedge is <j> = a - 0. This notation is chosen to be 
consistent with CH and the situation is illustrated in Figure 1 . 

We take ( x,y,z ) co-ordinates with x denoting the distance along the upper surface 
of the wedge, y the distance normal to the wedge face and z as the spanwise co-ordinate. 
We suppose that the corresponding velocities are (u,t), w) and denote quantities upstream 
of the shock and in the region between the shock and wedge by the subscripts u and s 
respectively. 

We assume that the fluid is a perfect gas so we have the equation of state 

P = ( 1 -l) C -Lpf , (2.1a) 

where p, p and T denote the pressure, density and temperature of the fluid respectively, 
whilst 7 is the ratio of specific heats and c p is the specific heat at constant pressure. Then 
the speed of sound a u in the upstream flow is defined by 

a u - (t 1 )^h> (2.16) 

Pu 
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where h u — c p T is the enthalpy of the fluid. It follows that the upstream Mach number is 
given by 

ft 

(2.1c) 


M. = ■+ 

d u 


The inviscid solution for hypersonic flow over a wedge is well known (see for example 
Hayes & Probstein ( 1966 )) for, if we define e = p u / ps, then we have 


7 — 1 

^= i2 — r ,(1 + 

7+1 


(7 — 1 )M 2 sin 2 <7 / ’ 


= 1 + 7M 2 ( 1 — e) sin 2 a, 

Pu 

= 1 + ~ 1 )( 1 - <?) M l si n 2 o-, 

flu Z 


( 2 . 2 a) 

( 2 . 2 b) 
(2-2 c) 

(2.2 d) 


tan <f> = e tan <7, 

( fi ll)„ = ^ coscr i («x)u = “C sin <7, - -eU sin <7, ( 2 . 2 c) 

where ti|| and u± are the velocity components parallel and perpendicular to the shock. 
Thus we see from ( 2 . 2 e) that the magnitude of the fluid velocity in the region between the 
shock and the wedge is given by 


U s = U cos <7 (l + e 2 tan 2 <7) 2 , 


( 2 . 3 ) 


which yields the slip velocity along the wedge. Thence, from (2.1c), (2.2c) and (2.3) the 
Mach number in the shock layer is given by 


M 2 = 


whence 


Ml 


2 M 2 cos 2 cr (l + e 2 tan 2 a) 

2 + (7 — 1)(1 — e 2 )M^ sin 2 cr’ 

2 Ml 


( 2 . 4 a) 


( 2 . 46 ) 


u 2 cos 2 a (l + e 2 tan 2 a) — (7 — 1)(1 — e 2 )M 2 sin 2 a 

Since we are interested in Tollmien-Schlichting instabilities of this basic flow we need 
to consider the governing equations which comprise the continuity, Navier-Stokes and 
energy equations for a viscous, compressible fluid. We non-dimensionalise these equations 
using the values of the flow quantities between the shock and the wedge. We write 


(x,y,z) = L(x,y,z), (u,v,w) = U 9 (u,v,w), 
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where the lengthscale L is the distance from the tip of the wedge to the position of interest. 
We non-dimensionalise time with respect to L/U s , pressure with respect to p s U ] and the 
other variables with respect to their values in the shock layer. If we define a Reynolds 
number for the flow by 

P s P 3 f 


Re = 


V* 


where p B is a typical viscosity then the flow is governed by the equations 

|+V.(^u) = 0, 


D u _ 1 

P~^I = -Vp+ 777 


and 


Dt 

DT 


Re 


2V.(^e) + V ((p' - |/<)V.u) , 


r2 D P 1 


dt ^ 7 Dt + PrRe " v ' ' ’ Re 

where p , p and T are the non-dimensional density, pressure and temperature. The shear 
and bulk viscosities p and p* have both been non-dimensionalised with respect to p s and 
Pr is the Prandtl number. Finally, the components of the rate of strain tensor are 


V.O VT) + (7 


(2.5) 

(2.6a) 

( 2 . 66 ) 

(2.6c) 


&ij 


du{ duj 
dxj dxi 


and the dissipation function $ is given by 

* = 2,.e : e + (p.' - (V.u) 2 . 

In the following analysis we shall assume that the Reynolds number is large. From 
(2.1a) the non-dimensional equation of state is 


lM]p = pT, 


(2.7a) 


and the non-dimensional enthalpy of the flow is given by 

h = T. (2.76) 

As in CH we shall restrict attention to the case when there is no heat transfer on the wedge 
surface so that dT/dy = 0 when y — 0 and the Prandtl is taken to be unity. Then the 
temperature of the fluid at the wall is given by 

Tu, = 1 + I (7 — 1) M 2 , (2.8) 
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where the subscript s has been dropped from the Mach number. 

The analysis presented in this paper is based on the assumption that the temperature 
dependence of the viscosity is described by Chapman’s law; i.e. 

fi = CT , (2.9) 

where C is a constant. The argument leading to the choice of scalings shown below was 
given by CH for the more realistic Sutherland’s law as well as for Chapman’s law. However, 
since the results presented in CH are restricted to the Chapman law scenario and since our 
attention is also focussed on this case we shall just give the scalings appropriate to this 
situation here. 

Then, following CH, we introduce the scalings 

a = SM* R e ~ ft and (7 - 1)MV = TS 2 , ( 2 . 10 a, b) 

where E and T are 0(1) constants. If (2.10a) is satisfied for a small then the shock will 
occur in the upper deck of the triple-deck structure and the effect of the the shock on the 
growth rate of the Tollmien-Schlichting waves can easily be determined. Consequently, we 
are assuming that <j > , the angle between the shock and the wedge face, is small. 

The problem is greatly simplified if the contributions from non-parallel effects can be 
neglected in a rational manner. This will be the case if the disturbance develops over a 
small streamwise distance. It can be deduced from the basic flow solution (2.2a) together 
with (2.4a) that e << 1 and for M \ > 0 in (2.46) we are forced to make the ‘Newtonian’ 
assumption 

(7 — 1) << 1, (2.11) 

if non-parallel effects are to be ignored. The scalings (2.10) take account of these conditions 
and these considerations furnish a lower bound for a. An upper bound for a is obtained 
by including the effects of nonlinearity in the lower deck which is the aim of the present 
study. In CH the linearised problem was obtained from the nonlinear problem they initially 
formulated. As in CH we intend to apply linearised shock conditions in the upper deck and 
this consideration whilst simultaneously demanding nonlinearity in the lower deck yields 
the desired upper bound on a. Specifically, the neglect of the non-parallelism and the 
imposition of nonlinearity imposes the constraints 

Re~~^ « a << Re~ (2.12a) 
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and so, from (2.10a) 

Re^° « M « Re&. (2.12 6) 

For a detailed description of the argument leading to (2.10) and (2.12) the reader is once 
more referred to CH. 

If the upper bound on a is satisfied then, as discussed in CH, we are able to ignore 
viscous effects in the entropy and shear waves which are inevitably produced when the 
pressure/acoustic wave meets the shock. With the scales described above and assuming 
that the shock has a viscous internal layer the entropy and shear waves will have y— scales 
much less than the width of the upper deck. 

From (2.106) if T w is given by (2.8) then 

T w ~ ^rE 2 <7" 2 >> 1, (2.13a) 

and to satisfy Af 2 > 0 in (2.46) we require 

TE 2 < 2. (2.136) 

The position of the shock is given by (2. 2d). Thus, in cartesian co-ordinates at the shock 
we have, from (2.2a) and (2.10a) that 

_ _3_ ^23 -j 

y = xe tan a m Re 16 M <*E x. (2-14) 

We re-emphasize the point here that much of the above is described in much greater detail 
by CH. 


§3. The triple— deck structure. 

It was shown by CH that the hypersonic flow described in the previous section is 
governed by a triple-deck structure. In deriving the relevant scales within each zone it 
is convenient to scale out a number of the variables. The triple deck structure relating 
to a hypersonic flow situation of the type we have here was deduced by Smith (1989) in 
the course of his investigations into the stability of supersonic flows. It is possible to use 
Smith’s results to infer the expansions which allow fi w , T w and M to be eliminated from the 
governing equations. These were given by CH but will be repeated here for convenience. 
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The x—, z- and t- scalings will be the same in each deck and are given by 


* = 1 + .Re-f CiriMt A-^X, 


(3.1a) 

z = Re~*C*T 2M~*Z, 


(3.16) 

and 

i = Re~* C*T W M? A~^r, 


(3.1c) 

where C = 1S the Chapman constant and A — Ax 2 , with A — 0.332, is 

boundary layer skin-friction from the undisturbed middle deck solution. 

the Blasius 

§ 3.1 The lower deck 



The lower deck is the region in which the nonlinearity of the problem manifests itself 
and the equations here necessarily contain viscous terms since it is this layer which is 
required to reduce the disturbance velocities to zero at the wedge surface. Indeed the 
relevant equations in the remaining zones of the triple— deck are linear as we shall describe 
presently. Within the lower deck we have scalings which take the forms 

y = Re-lciT$M*\-*Y, 


(3.2a) 

u ~ Re~*C*TwM* \*U, 


(3.26) 

v ~ Re-*C%TiM~$\W, 


(3.2c) 

w ~ Re-*C*TZM-*\*W, 


(3. 2d) 

p ~ 7 - 1 M" 2 + Re~* C* M~? P, 


(3.2e) 

T ~ T w , 


(3-2/) 

P~T~\ 


(3.2fir) 

We substitute these expressions into (2.6a, b) which give at leading order 


Ux + V Y + W z = 0, 


(3.3a) 

U T + U Ux + VUy + WUz = Uyy, 


(3.36) 

II 

o 


(3.3c) 
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Wr + UW X + V Wy + WWz = -P z + Wyy • (3.3d) 

The boundary conditions pertinent to system (3.3) require no-slip conditions on the wedge 
face so that 

U = V = W = 0 at Y = 0, (3.4a) 

and the solutions must match with those in the middle deck as Y — ► oo so that 

U — >Y + A(X,Z,t), W — >0, as Y — >oo, (3.45) 

where A is a displacement function which is to be determined. 


§5.j2 The middle deck 

The middle deck covers the extent of the undisturbed boundary layer. Since, from 
(2.13a), T w » 1 we require a thin transition region to reduce T to its free-stream value of 
unity. This transition region has thickness O ^(21ogT„,) - *^ in terms of the Dorodnitsyn- 
Howarth variable. Thus, the middle deck is itself composed of three regions: (i) a bound- 
ary layer region where T >> 1, (ii) the transition region and (in) the region above the 
transition region where T ~ 1. 

The thickness of region (i) is O ^Re~^C^T w S j and if we define 

y = Re-?C?T w y*, (3.5) 


then to match with the lower deck solutions as y* — > 0 we find that the solutions take 


the forms 

« ~ u* 0 (y*) + Re--*C*nM 1 *\--*Au* 0y ., (3.6a) 

v ~ —Re~* C* M - 2 A* A x Uq, (3.6 b) 

w ~ Re-*C*M~* A* (i?o“o) _1 D(X,Z,t), (3.6 c) 

p ~ 'y - * M~ ^ -j- Re~ * C*M~ * A * P, (3.6d) 

p ~ R* 0 (y*) + Re~*C*nM*\-*ARl r . (3.6c) 


Here Uq(?/*) an< i are undisturbed velocity and density profiles respectively and 

D x = — Pz- For Chapman’s law takes the form of a Blasius function for a compressible 
flow. 
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In region (ii) we write 

y = y c ~ Re-^C^T w (2\ogT w y^Y\ (3.7) 

where y c is a constant which determines the position of the transition layer. Here the 
solutions are very similar to (3.6) with the T w scales slightly altered for u and p. In 
particular, for u, Tj is replaced by ( 2 logT w )*. The solutions in region (in) are 

similar to those in region (i) with the simplification that Uq = Rq = 1 here. 

§ 3.3 The upper deck 

Finally, in the upper deck the basic flow quantities take their free-stream values. 
Thus, as in CH, the pressure/acoustic waves have scalings given by 


y = J&r*C*TjAr*A-*y, 

(3.8a) 

u ~ 1 + Re~*C* M~^ 

(3.86) 

( v,iv ) ~ Re~*C* M~* (u,u>) , 

(3.8c, d) 

p ~ 7 - 1 M -2 + Re-*C*M~% \*p, 

(3.8c) 

p ~ 1 + Re~*C<M^\^p, 

(3.8/) 

T ~ 1 + (7 — l)Re~*C*M%\^T. 

(3.8<jr) 


If we substitute (3.8) into (2.6) and (2.7) and rearrange the resulting equations we find 
that the pressure p satisfies 

PXX ~ Pyy ~ PZZ = 0. (3.9) 

The solutions in the upper deck must match with those in the middle deck as y > 0. 

From the y — momentum equation we obtain the boundary condition that 

p g = Axx at y = 0> (3.10) 

while p = P when y — 0. The remaining boundary condition is that 

p = 0 at y = y s , (3-11) 

where y 3 is the position of the shock. The analysis leading to this condition with our 
basic state is given in the Appendix of CH and will not be repeated here. The reader is 
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also referred to Moore (1954), Ribner (1954) and McKenzie & Westphal (1968) for the 
conditions which hold at the shock. This condition is obtained from the linear inviscid 
problem and will continue to hold for the nonlinear analysis presented in this paper since 
the nonlinearity is confined to the lower deck to the orders considered. 

For Chapman’s law the position of the shock is given by (2.14). Thus, from (3.8a), 
(2.13a) and (2.10a) we obtain 


Vs = 



At 

c is' 


(3.12) 


We shall now follow the convention of referring to the solutions of (3.9) as acoustic waves. 


§4. The nonlinear solution. 


In general the solution of equations (3.3), (3.4) and (3.9-11) is a fully numerical task 
but it is possible to make analytical progress by implementing a weakly nonlinear analysis 
of the system. The method used here is identical in spirit to that proposed by Smith 
(1979) who investigated the nonlinear stability of an incompressible Blasius boundary 
layer to Tollmien-Schlichting type disturbances. We differ slightly from Smith (1979) 
by considering a three-dimensional mode as in CH but this poses no formal difficulties 
whatsoever. 

Our objective is to monitor the streamwise development of a finite amplitude TS wave 
whose leading order wavenumbers in the streamwise and spanwise directions are a and /? 
respectively and whose frequency is Cl. Consequently, we seek modes whose fundamental 
component is proportional to 


E = exp 


+ (3Z - 



(4.1) 


If the weakly nonlinear disturbance is allowed to develop in the vicinity of the linear neutral 
point and if the relative amplitude of the disturbance in the lower deck is 0(h), h < 1, 
then following Smith (19(9) it is found that the amplitude A. of the mode will evolve on 
a O (h 2 ) lengthscale. However, there is a lower bound on the possible size of h which is 
derived by consideration of the non-parallelism of the flow. Using an argument formally 
identical to that of Hall & Smith (1984) it is straightforward to demonstrate that the 
following account is valid for all h in the range 0(Re~^ MtsT^) « h « 1. If the lower 
bound of this inequality is attained then the lengthscale over which the amplitude of the 
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perturbation modulates is identical to that over which the non-parallelism of the basic flow 
occurs. For full details of this aspect the reader is referred to Hall & Smith (1984). 

From (3.1a) we observe that the linear neutral stability point has been non-dimensional- 
ised so as to be at x = 1 and so we need to consider perturbations at the point 

x = 1 + h 2 x 2 - (4.2a) 

Since the skin friction A is a function of x it too will be slightly perturbed from its neutral 
value Ai according to 

A = Ai(l + h 2 \ 2 ). (4.26) 

Additionally, we choose to fix the frequency but allow the spanwise wavenumber to vary 
and write 

0 = 0i + h 2 0 2 . (4.2 c) 

These perturbations imply that for h << 1 we seek solutions of the lower deck equations 
(3.3-4) and the upper deck problem (3.9-11) in the forms 

3 

(U,V,W,P,A,p) = ((l + h 2 A 2 )F,0) + Y / h n (U n ,V n ,W n ,P n ,A n ,p n ) + 0(h 4 ). (4.3) 

n= 1 

Finally, the boundary condition (3.46) which describes the matching with the middle deck 
solutions now assumes the form 


U — > (l + h 2 X 2 )(Y + A), W — >0 as Y — ♦ oo. 


(4.4) 


To account for the slow modulation of the amplitude on streamwise lengthscales we 
introduce the co-ordinate X = h 2 X and then by use of multiple scales we replace all X 
derivatives throughout according to the recipe 


d 

dX 


d 

dX 


+ h 2 


,_d_ 

dX 


(4.5) 


Substitution of (4.3-5) in (3.3-4) leads to a hierarchy of problems at increasing orders in 
h. We address these problems in turn. 
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§^.l The first order problem 


The solution for the 0(h ) terms in (4.3) has the linear stability form 

(U u V u W i ,PuAuPi)= {u 1 {X,Y)^(X,Y),W l (X,Y),P l (X),A l (X),p l (X,y)) E+c.c„ 

( 4 . 6 ) 

where we recall that E is defined by (4.1) and here, as in the remainder of the paper, c.c. 
denotes complex conjugate. 

The problem for the unknowns in (4.6) is precisely the linear problem considered by 
CH. It is convenient to define 


£ = (*<*)* (l ) and £o = — » 

\ Oi ) (J 3 


and then we obtain 


aUi+faW! = — 

(ia)*A«'(£ o) 4, 


JJTTTT / Ai ^ d ^ 
3 -A* (Co) Jt o 


(4.7) 


(4.8) 


where Ai(£) denotes the Airy function. In order to satisfy the boundary conditions as 
£ > oo then 

=i5^jKPi, 


(4.9) 


with k = J^° Ai(q)dq. From the continuity equation (3.3a) we deduce that 


Vi 


(io:)3^4i 


K 


Ai'(0-Ai'Uo)-£ [ Ai{q)dq 
Ji o 


(4.10) 


Consideration of the upper deck problem (3.9-11) yields the solution 

-A\, for fil > a 2 , 


Pi 


a 2 sinh 

'(Pi 

i-a^Hvs-y) 

y/ i — ol 2 cosh 

\fil-a*)\y a 


(4.11a) 


and 


Pi = 


a 2 sin 

(a : 

2 -0l)Hys-y) 

y/ Ck 2 — /?! cos 

(a 2 - Pl)?y s 


A lf for fil<a 2 . 


(4.116) 


Note that the solution (4.116) for pi does not decay as y s — * oo. Since pi = P\ at y — 0 
we obtain from (4.9) and (4.11) the dispersion relation 

(/?i 2 - a 2 )* 


(ia) a /? 2 k 


^ z, (Co) tanh {{fi 2 — o 2 ) 2 y s ) 


for fil > a , 


(4.12a) 
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= for ft 2 <a 2 (4.126) 
S*'(to) tan {(a 2 — Si) 2 Ps) 

As mentioned by CH in the limit y s — ► 00 for j3j > a 2 the relation (4.12a) becomes that 
obtained by Smith (1989) for hypersonic flow with no shock. 

The neutral solution of (4.126) for a and 0 real is determined from the eigenrelation 

(/?i - at 2 )? = c 2 a? 01 tanh [{01 - <* 2 ) 3 y s ) , ( 4 -13) 

and £0 = — Cii'5 where Cj 2.3 and c 2 « 1.0. The corresponding expression for < a 2 
is similar to (4.13) and CH obtained neutral solutions to these two equations for various 
values of y s . The dispersion relation admits an infinity of solutions and in Figure (2) we 
present some of these solutions for the three cases y s = 1, 4 and 16. For each y s we 
observe that there is one mode for which a — ► 0 as 0 — > 00 and an infinity of modes 
all of which asymptote to a = 0 as 0 — * 00 . The Mach cone is defined by the line 
a = 0 so that above this line the acoustic waves are sinusoidal and supersonic whereas 
below this line the modes are subsonic. Hence for nearly all spanwise wavenumbers there 
are an infinity of supersonic modes and precisely one subsonic mode. We shall infer precise 
details concerning the asymptotic limit 0 — > 00 in Section 5 where we shall also discuss 
solutions of the amplitude evolution equation for A\{X) which is to be derived. 

§^.£. The second order problem. 

By examining 0(h 2 ) terms in (3.3) and (3.9) we conclude that the solutions of these 
equations at this order take the forms 

U 2 = U 21 E 2 + U 22 + U { 2 ? E~ 2 , (4.14) 

with similar expansions for V 2 , W 2 , P 2 , A 2 and p 2 . Here the terms U 2 \, U 22 , V 215 V 22 , W 2 \ 
and W 22 are functions of X and F, p 2 \ and p 22 are functions of X and y and the remaining 
terms are functions of X alone. 

From 0{h 2 E 2 ) terms in (3.3) we find, using manipulations similar to those described 
by Smith (1979), that 

otU 2 \-{-0W 2 i = — — f Ai) (ict) 3 / F(q)-\-Ai'{q) f Ai{t)dt dq-\-B 2 {X) f Ai{2?t)dt, 

\ / Jc o L Jto J 

(4.15a) 
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which, using the continuity equation (3.3a), leads to 


Vn = -2 [*(aU 2 i+0iW2i)dt . 
h o 


(4.156) 


In (4.15) we have made the definitions 


F(() = M2i {) [' Ai(t)R(t)dt, 

J 23£ 0 LAz(o) Joo 


<0 [^(y)]" 


where 


*(2*0 H -2-* [2i4i(0At M (0 + Az'(£ 0 )Ai'(0] • 

Imposing the boundary condition (3.46) as £ — > oo in (4.15a) gives 

•£ 

I *^£o 


(4.16a) 


(4.166) 


ia 2 / ~ \ 2 , y 00 r rS 1 r 00 . 

Aja^i = — rl 2 ^ 1 ) ( ia ) 3 / F + Ai'^) / Ai(t)dt d£ + B 2 / Ai(23<)df. (4.17) 

K K J J$ o •/£„ 


From 0(h 2 E 2 ) terms in the momentum equations (3.36, d) when evaluated at Y = 0 we 
obtain the relationship 


.. .2 ( d 2 U 21 , d 2 W 21 

M 3 + 


dt ; 2 


— 2ifi\P2i, 


(4.18) 




and elimination of P 2 between (4.15a), (4.17) and (4.18) yields a first connection between 
P 21 and A 21 . A second relation between these quantities is obtained by studying the upper 
deck problem (3.9-11) where, on considering 0(h 2 E 2 ) terms we find that 


d 2 p- 


21 


df 


- 4 (# - a )p 2 i = 0 


(4.19a) 


with associated boundary conditions 


dp 2 


dy 


P21 = 0 at y = y a , 


- = — 4o 2 A2i at y = 0. 


(4.196) 

(4.19c) 


Solution of (4.19), when combined with the matching requirement with the main deck 
solution that p 2 i = P 21 at y = 0 yields the relation 


P 21 = 


2a 2 A 2 i 


r tanh ( 2 (/5i - a 2 ) » y 3 ) , 


(4.20) 
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for {d 2 > a 2 with an analogous expression when (d 2 < a 2 . This second expression for 
P 2 \t A 2 \ enables B 2 (X ) to be found, as previously described. The resulting expressions 
are complicated and so we shall not state them here. 

Turning to the 0(h 2 ) mean flow terms U 22 , V 22 , W 22 , P 22 , A 22 and p 22 in the governing 
equations we obtain the results that 

V 22 = 0, (4.21a) 



This latter equation has solution 


(4.216) 


aU 22 + fii fb* 22 


(*«) 


_ 2 9 
3ft 


Ai 


K 



f**(q)dq, 


(4.22 a) 


where 


/"«) = >• i (Ai(0) M 


Ai'(0~Ai'(& 0) — € / Ai(t)dt 

0 


+ C.C. 


(4.226) 


Applying the boundary condition (3.46) as £ — » 00 implies that 

a(m)-t,ii '°° r* 


A22 — 


! £bl|— p f rim. 

M * Jz 0 J 00 


(4.23) 


Finally, the upper deck problem (3.9-11) yields p 22 = P 2 2 = 0. 


§^.5 The third order problem. 

At this order we determine the amplitude equation for the unknown function A\(X). 
This equation arises from a solvability condition on the 0(h 3 E) terms in the systems (3.3) 
and (3.9-11), solutions of which are sought of the form 

U 3 — EU 31 + E 2 L t 3 2 + E 3 Uzz + 6^34 + E x Tj£ + E 2 Uw + E 3 U^ 3 \ 


with similar expansions for V3 , W3 , A 3 , P 3 and p 3 . The desired solvability criterion is 
obtained from consideration of the coefficients of 0 (h z E ) terms in the lower deck equations 
(3.3) and this leads to 

^-(at/31 + PiW 31 ) - +0 iW 31 ) = + (4.24 a ) 
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where 


O 

-G((,X) =(to)-*({ - (o Waft + A Jfc) + (ia)*(£ - f 0 )Aj(ot>, + AW,) 

uJ\ 

+ a\2V\ + iflifaPi + 1 + /?i Wi)^(c*{72i + P 1 W 21 ) 


o 

+ i(aU22 + A W22)(<*J7i -f ftWi) + (ia)*Vi ^ ( 0^22 + /?i W 22 ) 

+ M^/ C) ^(a^2i +A^2i) + ((ia)i) (C) V 21 ^(a^ (c) + /W (c) ). 


(4.246) 


The boundary conditions obtained from (3.4) are 


U 31 (Zo,X) = W 3 i(fo,X) = 0, 


(4.25a) 


U 3 i(oo,X ) = A 31 (X) + A 2 A 1 (X), W 31 (oo,X) = 0 . 


(4.256) 


By evaluating the streamwise and spanwise momentum equations (3.36, d) at £ = £ 0 we 
get 


i*a~*PlP 3 1 + ita-SfafoPi 


. 1 _ 2 


-( 


d 2 U 3 1 d 2 W 31 

« + /?i 


) 


(4.26) 




d£i ' ^ d£ 2 

The homogeneous forms of equation (4.24) and boundary conditions (4.25) are sat- 
isfied by aUi + d 1 14 7 i and so (4.24-5) only has a solution for a- U 3 1 + 0\W 3 \ if a certain 
compatibility condition holds (as in Stuart (1960) and Watson (I960)). This condition is 
derived by considering the adjoint system of the problem as in Hall & Smith (1982). We 
multiply (4.24a) by the function 


a'«) - *«) - gjMcia 


(4.27a) 


where 




Uo \Ai{q)] 2 

and integrate over the range ((0, 00). This gives 


(4.276) 


i^a-^AiUo)^P 3 i - 

L (Co) 


^^aA 2 ij + (.«)“* J” K’(t)G(t,X)dt + J” K{$ d£ 


(4.28) 
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In the upper deck the problem for p 3 i takes the form (on equating terms of 0(h 3 E) 
in (3.9)), 

^ _ (f _ a 2 )P3i = 2io§ +2A&P,, (4.29a) 

oy * oX 

with boundary conditions 

P31 = 0 at y = y a , (4.296) 

^4 = -q 2 A 3 i + 2 iot^- and p 3 i = P 31 at y = 0. (4.29c) 

dy 8X 

Solution of (4.29) yields a second expression for 

i* a~* Ai(£o)0i P 31 + 77 jpr Q d. 3 i 

Uoj 

in terms of the solutions derived when considering the first and second order problems 
in Sections 4.1 and 4.2. Then, after considerable manipulation, substitution of (4.246) 
for G(£,X) in (4.28) enables the desired evolution equation to be obtained. Since the 
expressions occurring in this equation are long and complicated we have relegated these 
details to the Appendix. In the following section we analyse some solutions of the amplitude 
equation and comment upon some of the implications. 


§ 5 . Results and discussion. 

By implementing the procedures described in the previous section we obtain a classical 
evolution equation for the amplitude Ai(A’) which takes the general form 


— = (ui/?2 + a2^2)Ai + a$A\ |.4i | 2 , (5-1) 

dX 

where the definitions of the complex valued constants d \ , 0,2 and ct 3 are given in the Ap- 
pendix. As is usual in this type of weakly nonlinear analysis interest is primarily centred 
upon the sign of the real part of a 3 . It is easily deduced from (5.1) that 

4 - (li^ 2 ) = [2Re(a 2 \ 2 +a 1 f3 2 )]\A 1 \ 2 +2Re(a 3 )\A 1 \ 4 , (5.2) 

and thence if Re(a 3 ) > 0 nonlinear effects are destabilising and the equilibrium state with 
amplitude i 

\ -fe(°3 I P (5 .3) 


Re(a i/?2 + 02^2)] 
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is unstable. Conversely, for Re(a 3 ) < 0 nonlinearity tends to stabilise a linearly unstable 
mode. (Here i?e($) denotes the real part of the complex number 4>.) 

We now illustrate solution characteristics for the values of y s considered previously. 
Firstly, in Figure (3) we show the dependence of the real part of the coefficient a j as 
functions of (3\ and where a is then given by the solution of the eigenrelation (4.12). We 
see from these diagrams that the behaviour of ai is markedly different for the subsonic 
mode as opposed to the remaining modes. Specifically, for sufficiently large y s the real 
part of cii appears to be negative for all and for all modes except the subsonic one. 
Figure (4) shows that the constant a 2 has the property that Re(a ,2 ) < 0 for all values of y s 
and /?i. In particular, the subsonic mode has the greatest absolute value for fje(a 2 ) and 
the ratio between the size of this constant for the subsonic mode to that for the supersonic 
modes increases as the shock position y 3 moves away from the wedge surface. 

When the amplitude equation (5.1) is linearised we observe that the sign of Re(a 2 A 2 + 
ai/? 2 ) determines whether small amplitude Tollmien-Schlichting waves are stable or un- 
stable according to a usual linear stability analysis. As previously described, it is the sign 
of Re(a 3 ) which governs the effect of the nonlinear term on the stability of near neutral 
modes. Crucially, Figure (5) demonstrates that for all y s and /?i we have Re(a 3 ) < 0 and 
hence for each mode the effect of nonlinearity is stabilising. The size of the nonlinear term 
decreases as y s increases and for sufficiently large y s the nonlinear term relating to the 
subsonic mode far exceeds in value those for the remaining supersonic cases. 

Above we have provided a brief description of the results of the weakly nonlinear 
calculations for a selection of shock positions y s . We now turn to consider some analytical 
results which are derivable in certain important limits. CH identified three particular cases 
in which asymptotic progress is easily accomplished. Firstly, we consider the case of a weak 
three-dimensionality effect in the TS disturbance (4.1), i.e. f3 3 << 1. Then the solution 
of eigenrelation (4.12) takes the form 


a ~ 


(n + |)7T 


Vs 


+ 


y» 


2(n + |)7r 



(5.4a) 
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where n = 1,2,... and where we recall that C 2 ~ 1.0. Then, using the functions defined in 
the Appendix it follows that at leading orders in 0i the evolution equation for A\ becomes 



where the constants Lj are defined in the Appendix. On evaluation this amplitude equation 
becomes 


dAi 

dX 


A 20\ 


,yl(n+ \)i 


-(-0.378,0.273) + 

Pi 


P1P2 


y»{n + *)t 


(-0.001,-0.932) + 


ys 


(n + |)3 


r (0, 0.318) 


A\ 


+ —(—0.279, 7.03 x 10 -2 )Ai|Ai| 

Vs 


(5.46) 

We may deduce from this equation that for 0i << 1 all the coefficients in this amplitude 
equation become small. Further, on considering the coefficient of the linear term we see 
that variations in the perturbed spanwise wavenumber 0 2 have little effect on the linear 
growth rate of the TS wave. Finally, for 0 2 = 0 we have the possibility of a finite amplitude 
stable mode of size 

l.lQyf 


Ac = 


)ys AT 
ttv-a2 


(n + 2)3 

A second case for investigation is that of supersonic modes with a » 1 . From Figure 
( 2 ) we observed that such neutral modes have "w a and, more accurately, we find that 


a = Pi + I — ) 


1 


y a J 2 /?! 


-i+Y 


+ 


A(-) + 


(5.5a) 


ys ) S0f ' C 2 y s \ys J 

where n = 1,2, Thus for 0\ » 1 the possible a values for neutrally stable linear TS 

waves are separated by 0(l/n). In this case the amplitude equation asymptotes to 


PjylMio) dAi _ 
(nn) 2 dX 


\ 2 0i L 2 + 02 


( PML±\ 
V ( n7r ) 2 J 


: 0?(c 2 L 5 + 4 L 6 ) : , 1 , 

- 4l + -+ 2 I 7 +4L s ) 


(5.56) 
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Performing the necessary evaluations yields 


dAi 

dX 


(-8.016,5.781) 


n 


% 


10 

3 


y» 


-A2 + i/3 2 


Ai +(-6.058,-0.377) 


y 3 . 


-Ailiil 2 , (5.5c) 


and hence we can obtain supercritically stable modes of amplitude « 1.15 /? x *\/— ^2- We 
note that in this a, /5 1 limit the scaled equilibrium amplitude is small and is independent 

of the perturbed spanwise wavenumber /?2- 

A final regime for asymptotic investigation is that of the solitary subsonic mode in 
Figures (2), i.e. those modes for which a — > 0 as /?i — * 00. The solution of eigenrelation 
(4.12a) in this case assumes the form 


Q ~ 


1 

I^M 3 


+ 


(5.6a) 


and then the amplitude equation (5.1) simplifies, at leading order in 1 to 


C 2 L 1 — 2^4*(£o) dA\ 
W4 ~dX 


L2 \ ^2(^2 713 + L4) 

WJ 


M + 


(c 2 L 5 + 2 I 6 ) x 1 a 12 

0^(0217 +2L 8 )^ l! 11 ’ 


or, 


(1 A 1 

-4 = [(-0.376, 0.117)A 2/ 3- 3 +(0.125, -0.373)^ 2 ^r 4 ] A x 

dX 

+ (-0.321,-8.83 x 10” 2 )i8 1 _5 Ai|i 1 | 2 . 


(5.66) 


This evolution equation admits a stable finite equilibrium solution in which A\ = 0(/3 1). 

We remark that the three asymptotic limits briefly alluded to above agree well with the 

numerical solutions sketched in Figures (3)-(5). It is straightforward to make a comparison 

of the relative likelihoods of the equilibrium modes in the (3 1 >> 1 case. We have seen that 

_ 1 

for f3i >> 1 the supersonic modes have amplitudes 0(^ 1 3 ) whereas the corresponding 
amplitude for the subsonic mode is much larger, 0(/3 \ ). Hence, in practice, at large 
spanwise wavenumbers /?i the subsonic mode will likely dominate the other modes and is 
therefore much more likely to be observable. 

To conclude, we have presented an account of some properties of classical weakly 
nonlinear stability theory as applied to an, admittedly idealised, hypersonic flow problem 
originally considered by CH. Our main conclusion is that for all the parameter regimes 
investigated near linear neutrally stable modes appear to be supercritically stable; i.e. the 
nonlinear effects tend to stabilise the slightly unstable linearised modes. 
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There is much further work which is possible in relation to this problem. The most 
obvious is that of a full nonlinear solutions of the lower deck equations (3.3) together 
with boundary conditions (3.4); work which when combined with the relevant upper deck 
problem clearly needs to be numerical in character. A second course of action is that of 
considering viscosity/ temperature dependences other than that of Chapman’s law. As 
mentioned earlier CH and Blackaby et. al. (1990) consider situations in which the more 
physically realistic Sutherland’s law is used although CH did not calculate any results for 
this case. The effect of changing the viscosity law on our results would be of great interest 
for, of course, if it transpires that the properties of our solutions are quite sensitive to the 
actual laws invoked then this raises the question as to which to use in future studies. 

One complete subject area which has been deliberately excluded in both CH and the 
present study is that as to the nature of any nonlinearities which might occur in the vicinity 
of the shock layer position y s . The scalings chosen here were such that (2.12a) holds: 
technically once a reaches a value of 0(Re~^») then the shock layer problem becomes 
nonlinear itself for we can no longer ignore viscous effects in entropy and shear waves 
which are produced when the acoustic wave meets the shock. Finally, w^e mention the 
observation of CH who concluded that in addition to the neutral modes sketched in Figure 
(2) there are other modes with zero growth rate but with infinite frequency and therefore 
these modes cannot be described within the asymptotic framework used here. Investigation 
of both the linearised and nonlinear problems which arise within the necessarily modified 
structure would be extremely worthwhile. 

There are, therefore, many directions in which the results presented here may be 
extended. Progress in any of these would be valuable indeed. 
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APPENDIX 


In this appendix we summarise the results of manipulating (4.28) and using the first 
and second order solutions in order to derive the desired amplitude equation. We recall 
that aU\ + pi W i is given by (4.8) and by considering the 0(h) terms in the lower deck 
equation (3.3 d) we obtain W\. Similarly, W 21 can be found by examining 0(h 2 E 2 ) terms 
in this equation and then U 21 derived via (4.15a). Finally, V\ and V 21 are given by (4.10) 
and (4.156) respectively and thus <?(£, X) deduced from its definition (4.246). We may now 
proceed to substitute in (4.28) and the resulting classical Stuart- Watson type evolution 
equation for Ai(X) is as follows: 

If p 2 > a 2 then 


<«*£. - ( (2ft a2 a ) tanh( x //??-^y.) 

(Pi — a 2 ) 2 \ at* V 


a 


1 ( L$ + 
+ a 3 ' 


L 7 + $>+Lg 


2 _ „,2' 

a 2 

2 _Pb 

(Pl-cc 2 )* 

iiliil 2 , 


y/P\ - « 2 


Vs 


cosh 2 ( y /pf^a 2 y s ) 


dA\ 

~dX 


*L 2 A 2 + ^£3 + , ° ^ , Li (tanh( i/ft-a'y.) 


\/fi{ - a 2 Vs 


cosh 2 ( \fpp~^ a 2 y s ) 


A\ 


( A.la ) 


where 


_ 2a 3 /3 2 




v^F - « 2 


t anh( 2 yfpf- « 2 y s ) , 


and where are complex constants defined presently. If however /? 2 < a 2 then 


d L, + _ V EZEv± 

(a 2 - P\)l \ a * cos 2 ( v / or 2 - /?Fj/ a ) 


dAi 

~dX 


5 a 3 Qo 

a*L 2 \ 2 + —^L 


+ OL 


Pi (a 2 — 

7 ( L$ + 

1*7 + ^ — -^8 


a 2 - Pi)* \ ^ cos 2 ( v /a 2 - t/a) / 


iiliil 2 , 


(A.16) 


where 


= / 2 °2 =^f tan ( 2 \/ a2 ~ £?&)• 
V 1 “ - /h 
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The constants here are defined by 


L x 


(Tl) 


iAi'Uo) 


K(OM() 



(A.2a) 


(Tl) 2 iiAi((o)Ai'(( 0 ) 


Li = — 


K 


L, = 


izAi'do) 


At(6) 


K. 

A ds 


m)Ai (o ( f -A-i f Ai(t)dt ] d( 

'to \Jt o [ 4 t(s)J -/«> / 


Z74 — 

(T20)(T14) + (T21)(T7) 


K\K\ 


(A- 2b) 

(,4.2c) 

(A.2d) 
(A. 2e) 


and 


r (T20)(T13) - (T21)(T12) 

— 1 .19 



(-4-2/) 


(4.2-7) 


(A.2h) 


The complex constants Tl, T 7, T12, T13, T14, T20 and T21 are defined and evaluated 
in the Appendix of Bassom (1989). All the above expressions were found numerically. 
Initial conditions at a suitably large value of £, say were found using asymptotic 
methods and the defining equations integrated using a fourth order Runge-Kutta scheme. 
The integrations necessary in the evaluation of (A2) were performed using the trapezium 
rule combined with Richardson extrapolation. Various values of £oo and step length were 
used and the results quoted below are believed to be accurate to within 0.1%. Further 
checks were made on these computed values by using the program to calculate some of the 
coefficients given by Hall &; Smith (1982) and defined in their equation (2.27). Using the 
above techniques we found that 


Ti « (0.4766,0.9858), T 2 « (0.0277,-1.202), 
T 3 « (1.916, 1.450), T 4 « (-0.9573. -0.7247), 
T 5 « (-0.8346, 1.389), I 6 « (2.092. -0.9825), 
L 7 « (-4.633, -0.6513), I 8 « (0.3473, 1.536). 
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Figure (1). The geometry of the wedge and shock for a high Mach number flow. 





Pi 


Figure (2). (a) The first three neutral curves a = a(/3 1 ) given by (4.12) for y 3 = 1. (b) 
The first seven neutral curves for y 3 = 4. fcj The first ten neutral curves for = 16. 
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Figure (3). The real part of the coefficient as defined in (5.1). (a) The values 

taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the seven 
curves in Figure (26) where y 3 = 4. (c) The values taken on the ten curves in Figure (2c) 



Figure (3). The real part of the coefficient ) as defined in (5.1). (a) The values 

taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the seven 
curves in Figure (26) where y s = 4. (c) The values taken on the ten curves in Figure (2c) 
= 16. 


Figure (3). The real part of the coefficient ai(/?i) as defined in (5.1). (a) The values 
taken on the three curves in Figure (2a) where y 3 — 1. (b) The values taken on the seven 
curves in Figure (2b) where y 3 = 4. (c) The values taken on the ten curves in Figure (2c) 
= 16. 
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Figure (4). The real part of the coefficient a 2 (/?i) as defined in (5.1). (a) The values 
taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the seven 
curves in Figure (2b) where y s = 4. (c) The values taken on the ten curves in Figure (2c) 
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Figure (4). The real part of the coefficient a 2 (Pi) as defined in (5.1). (a) The values 
taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the seven 
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Figure (4). The real part of the coefficient 02 (/?i ) as defined in (5.1). (a) The values 
taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the seven 
curves in Figure (2b) where y s — 4. (c) The values taken on the ten curves in Figure (2c) 



Figure (5). The real part of the nonlinear coefficient as defined in (5.1). (a) The 

values taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the 
seven curves in Figure (2b) where y s = 4. (c) The values taken on the ten curves in Figure 
(2c) where y s = 16. 



Figure (5). The real part of the nonlinear coefficient a 3 (/3\) as defined in (5.1). (a) The 
values taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the 
seven curves in Figure (26) where y 3 = 4. (c) The values taken on the ten curves in Figure 
(2c) where y s = 16. 



Figure (5). The real part of the nonlinear coefficient az(fli) as defined in (5.1). (a) The 
values taken on the three curves in Figure (2a) where y s = 1. (b) The values taken on the 
seven curves in Figure (2b) where y s = 4. ( c ) The values taken on the ten curves in Figure 
(2c) where y s — 16. 
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